Statistics of skyrmions in quantum Hall systems 



Abstract 

We analyze statistical interactions of skyrmions in the quantum Hall system near a critical filling fraction in the 
framework of the Ginzburg-Landau model. The phase picked up by the wave- function during an exchange of two skyrmions 
close to v = f/(2n + f) is n[S + l/2(2rt + 1)], where S is the skyrmion's spin. In the same setting an exchange of two 
fully polarized vortices gives rise to the phase n/(2n + 1). Skyrmions with odd and even numbers of reversed spins have 
different quantum statistics. Condensation of skyrmions with an even number of reversed spins leads to filling fractions 
with odd denominators, while condensation of those with an odd number of reversed spins gives rise to filling fractions 
with even denominators. 
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^ , Recently skyrmion excitations near the filling factor v — 1 have drawn remarkable attention jl], |^] . There are two kinds 
topological excitations in single-layer quantum Hall systems. When the system is fully polarized, the relevant charged 

G|uasiparticles are topological vortices. At v = l/(2n + 1) the charge of such a vortex is e/ (2n + 1) and the spin is 1/2 (2n + 1). 
has been suggested Q that for a weak Zeeman coupling the lowest energy quasiparticle is a skyrmion or a slowly varying 

topological texture. Its electric charge is still e/(2n + 1) but its total spin S can be substantially larger. Such large spin but 

^moderate charge excitations can explain the observed depolarization as the filling factor slightly deviates from v = 1/ (2n + 1 ) 

■ It has been proposed to describe skyrmion excitations in terms of an effective ferromagnetic model |1[ . The model explains 
C^he dependence of the skyrmion size on relative strenght of Coulomb and Zeeman interaction. However, it does not predict 
correctly skyrmion energies. The energy of skyrmion and antiskyrmion seem to be the same. Microscopic calculations show 
that it is not the case. The effective model has been derived jl], ^ from the Ginzburg-Landau model of the quantum Hall 
system |], In the derivation, among other things, the charge density was assumed to be constant. This leaves polarized 
vortices outside the reach of the effective model. In terms of the Ginzburg-Landau model the difference between a vortex 
and a skyrmion is that a vortex is fully polarized while inside the skyrmion core there is some number of electrons with 
reversed spins. Such a depolarization helps to minimize the total Coulomb energy. It costs some Zeeman energy but it is not 
an obstacle when the Zeeman coupling is sufficiently weak. 

In this paper we consider some topological properties of skyrmions and vortices in the framework of the Ginzburg-Landau 
model. We show that the Magnus force acting on vortices and skyrmions is the same. We also consider statistical interactions 
of skyrmions. We will show that during a clockwise exchange of two skyrmions the wave-function picks up the geometrical 
phase 

r =t s+ whr) ] l'' s- 4 * <"-««'■ (1 » 

where S is the skyrmion's spin. This formula applies also to vortices when we put S — S vortex = (2n+i) 1 Thus the Ginzburg- 
Landau treatement shows that the statistical interaction is much stronger for skyrmions than for vortices. It is proportional 
to the number of reversed spins. 
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1 The Ginzburg-Landau model and topological defects 



We wish to investigate the statistical interaction of skyrmions in the quantum Hall effect. We will describe this system by 
the Zhang-Hansson-Kivelson model [^) generalized by Kane and Lee JtJ to describe unpolarized quantum Hall systems. 
The Lagrangian of the model is 

L = i^(d t + ia )<P-^ | (dk + iiak + eAkj^l 2 + ^ e ^^« CT ~ ) -p ] 2 +7[0!01 - 0^t] +m[0^1 + ^^t] • (2) 

Greek indices denote 0,1,2, while Latin indices take values 1,2. When they are repeated, summation is understood. We 
use the signature (+,—,—). (f> is a two-component complex scalar field, <fi = a M is a statistical gauge field while 

Ak = —\Bskix l is a gauge potential of the external uniform magnetic field B directed down the z-axis. epo is the positive 
background charge density, related to the external magnetic field by eB = 20po- For the boson field to represent fermions 
the parameter O must take one of the values = (2n + l)7r, where n is a nonnegative integer, m is the effective electronic 
mass and 7 is the effective Zeeman coupling. /1 is a chemical potential, chosen so that the ground state of the system (|J) 
is, up to a gauge transformation, the solution 4>i — i/po,4>-\ — with the statistical gauge field "screening" the external 
magnetic field, ak + eAk = 0. The chemical potential has to be ji = —7. In this ground state the system is fully polarized. 
The Lagrangian (Q) should be supplemented by the Coulomb interaction term 

-\f d 2 x'[^m - po] f /s frt^-/) _ Po ] , ( 3 ) 

* J I x — x I 

where e is a dielectric constant of the host material. 

There are two types of relevant topological excitations in the model. One of them is simply a vortex or a fully polarized 
quasihole. This configuration is given by the Ansatz 

01 = U(r)e- ie , 
eB 

a g = —r + a v {r) , 

a = b v (r) . (4) 

The modulus f v (r) has to interpolate between f v (0) = and f v (oo) = i/po- The asymptote of a v (r) at infinity is a v (r) « -. 
The vortex is by definition restricted to the lower spin component. This restriction is consistent with field equations of 
the model (^) for any value of the Zeeman coupling 7. Such a solution indeed exists as is well known from the studies on 
vortices in the Ginzburg-Landau model of the fully polarized quantum Hall effect j^]. The total charge is proportional to 
the statistical magnetic flux and equals to e/(2n +1). The distribution of the electronic density is such that it vanishes in 
the very center of the vortex and tends to — epo ouside the vortex core. This is not the most favourable charge distribution 
from the point of view of both the Coulomb interaction and the quartic self-interaction term. While the total charge is fixed, 
the Coulomb interaction tends to minimize the integral of the charge density squared. In other words it tends to make the 
charge distribution as diluted as possible. When we restrict to the spin-down component, the deviation of the charge density 
from — epo must equal to epo at the center of the vortex. This is a topological obstruction. At the same time the total charge 
is quantized so there is not much freedom left to make the charge distribution more dilute. However, for sufficiently small 
Zeeman coupling the configuration may find a way to minimize its Coulomb energy by nucleation of the spin- up component. 
A more general configuration, with more freedom to minimize the Coulomb energy, is 

h = g{r) , 

eB 

ao = —r + a s [r) , 

a e = b s (r) . (5) 

Now we admit a nonzero spin-up component. We want it to be nonzero in the vortex core, 3(0) ^ 0, and to cost a minimal 
gradient energy so the phase of the spin-up component has to be constant. For a finite energy configuration 3(00) = 0. 
Excitation of the upper component costs some Zeeman energy so it is not energetically favourable except for small effective 
Zeeman couplings 7, where the gain in Coulomb energy can prevail the loss in Zeeman energy. Note that with the Ansatz 
(|^) the average spin points up in the middle of the soliton while it points down outside the core. This configuration is a 
skyrmion but in the language of the untruncated model (||). As compared to a uniform background with charge density po, 
the | component adds a negative charge Qf, 

/•OO 

Q T = -27re / rdr g 2 (r) , (6) 
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while the | component contributes a positive electronic charge deficit 

Q l = -2ne / rdr [/ 2 (r) - p ] . (7) 
Jo 

The total charge is the same as for the fully polarized vortex; the two contributions add up to + Q i = ( 2 n+i) • 

Thus there is an interplay of two factors, namely the Coulomb and the Zeeman energy. For strong Zccman coupling 
vortices are the relevant quasiparticles. At weak Zeeman coupling vortices are still solutions of field equations but they have 
higher energy than skyrmions. In this limit skyrmions become the relevant quasiparticles. 



2 Magnus force acting on skyrmions and vortices 

We apply the adiabatic approximation to investigate the dynamics of vortices and skyrmions. The only terms in the La- 
grangian (j^) , which can contribute to the terms in the effective mechanical Lagrangian which are linear in velocity, are 

I%) f = J d 2 x [i$d t ct> - ^s kl a k d tai } . (8) 

The prescription for the adiabatic approximation is as follows. Take the static vortex or skyrmion solution described by the 
fields {4>(x — £), a^{x — £)} and located at an arbitrary position £. Promote the parameter £ to the role of a time-dependent 
collective coordinate In this way the static fields become time-dependent. The final step is to substitute such a time- 
dependent field configurations to the Lagrangian and integrate out their spatial dependence. After the spatial integration 
one should be left with a purely mechanical Lagrangian being a functional of the trajectory 

The first term in does not contribute to the effective Lagrangian. The static field is in the Coulomb gauge, d k ak = 0, 
so that the gauge field can be expressed as a k = EkidiU, where U is an auxiliary potential. In the adiabatic approximation 
we replace U(x) — + U[x — The second term in (g) becomes 

/ ^^ £kl£ ^ d ^ Ud ^ £ '^ U ^ = J d2x [-^mnd m Ud t d n U] . (9) 

It can be shown, by an integration by parts, that this term does not contribute to the effective Lagrangian. Thus the only 
contribution to the part of the effective mechanical Lagrangian which is linear in velocity comes from 

L eff = J d 2 xic^d t (j) . (10) 

When we introduce moduli and phases of the scalar fields, 4>A — ^/PA~e lXA with A =|,i, the effective Lagrangian will split 
into contributions from spin-up and spin-down components 

L ( Sf = - J d 2 x [p T 9 tXT + PldtXl] • (11) 

Let us apply the procedure to a single vortex or skyrmion. The static configuration is described by (|^). In the vortex case 
we can put formally g(r) — 0. Let us take the trajectory £(t) which passes through the origin at t — 0. For a very small £ 
the fields in ([ll]) can be expanded as 

PA {x-i) = PA (r) - (e cos9 + e smfl)^(r) + 0{\ £ | 2 ) , 

d t XA(x-0 = -n A (-e^+e C -^) + O(\Z\) ■ (12) 
With this expansion the effective Lagrangian becomes 

L« = -ire ki ei l E "a\Pa{°°) ~ Pa{0)] + 0(\ £ | 2 ) . (13) 

A 

The system with just one quasiparticle is translationally invariant. The term proportional to £ki£, k £, 1 is the most general term 
linear in velocity which is, up to a total time derivative, translationally invariant. This proves that the term 0(| £ | 2 ) in 
Eq. (|l3|) vanishes identically. 

There are two contributions to the effective Lagrangian. The contribution from the spin-up component vanishes because 
P]{oo) = ri|P|(0) = 0. Thus the only contribution comes from the lower component ( pi(0) = , pj(oo) = po ) and amounts 
to 

L^=^P mei l ^^ T) e kl ^ 1 . (14) 
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The Magnus force is exactly the same for skyrmions and for vortices. It equals the Lorenz force acting on a particle with the 
charge 2 n+i • ^ s discussed in || this is not a mere coincidence. A quasiparticle is not an independent object but a defect 
composed out of electrons. As the location of defect moves it act with a Magnus force on the surrounding electrons. They in 
turn respond with a current which interacts through the Lorentz force with the external magnetic field. Thus the two forces 
combine to be just one. 

If the inertial mass were zero then the Magnus force would prevent a quasiparticle from moving with respect to the 
condesate. Note that in our derivation we have been slowly varying the location of a quasiparticle without moving the 
condesate so that the £'s in Eq.([l4j) are quasiparticle coordinates with respect to the condesate's frame. The Lorentz force 
acts on the condensate as a whole. Namely, if the fields { (f>{t, x), a^(t, x) } are solutions of field equations of the model (^), 
then there are also the following boosted solutions 

cf>(t, x) = <t>[t,x- R{t)]e lXB , 

a a (t,x) = a [t,x-R(t)] - R k (t)a k [t,x ~ R(t)] , 

a k (t,x) = a k [t,x- R(t)} , 

Xb = mR k x k - -m[R k R k ]t - e / * dr R k (T)A k [x - R{t)} , 
2 J t 

R(t ) = , 
m 

A k {x) = -—e H x l . (15) 

Any solution can be boosted to move along an electronic cyclotron orbit. Even the uniform condesate feels the Lorentz 
force in spite of the fact that the external magnetic field seems to be screened by the statistical gauge field. If there is a 
quasiparticle in the original solution then after the boost it follows the cyclotron orbit but it does not move with respect to 
the condensate. The Magnus force remains zero. 



3 Statistical interactions of quasiparticles 

We have shown that the Magnus force is the same for both skyrmions and vortices. Now we proceed to their mutual statistical 
interactions. We will show that the statistical interaction of skyrmions depends on the total number of inversed spins they 
carry. 

Once again we will apply the formula (|ll|) but this time to the system of two distant antiskyrmions. The total phase of 
the lower component reads 

Xi {x) = -Argix-b) - Arg(x - &) (16) 

while the upper component phase is constant. Thus, once again, only the J, component contributes to the effective Lagrangian. 
The formula ( p"l| ) simplifies to 

L eff = ~\ SX P/^Xi] > 

where Spi = p[ — po is a deviation of the down component's density from the uniform background. When the distance 
between the skyrmions is large as compared to their widths the deviation can be approximated by the sum 

S Pl (t, x ) w S Pl [x - 6 (*)] + 6 PL [x - £ (t)} (18) 

of two nonoverlapping rotationally-symmetric deviations. When the formulas (|l6| ) and ( |l8| ) are substituted to (|l7|), there 
appear four terms to be integrated out. Two of them are " self- interaction" terms which lead to Magnus forces as discussed 
in the previous section but in addition there are two mutual-interaction terms 

l£} f = J d 2 x {S Pl [x - 1*1 (t)]d t Arg[x - £(*)] + 6 Pl [x - i 2 (t)]d t Arg[x - ^(t)}} . (19) 

For very distant skyrmions the density deviations can be approximated by 5f>i[x ~ £i(i)] ~ —^5^ 2 \x — where Qi is 

the charge deficit of the spin-down component, compare with (^). In this approximation the integral ( Jigj ) becomes 

L eff = "T / ^ m ^- UmtArg[x- £(*)] + 5^[x-Ut)]dtArg[x-Um = -^^Argfo® - &(*)] , (20) 
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which is just the statistical interaction we were looking for. The prefactor Qi/e is the total number of electrons missing in the 
lower component as compared to the uniform ground state. As the number of electrons in the upper component is — Qj/e, 
the total spin of the skyrmion as compared to the uniform background is S — \ {Q] — Qi)/e. Qj = — Qi — 2 n+i > so the 
statistical interaction becomes 

^ = -[^^^i^ 9[m - m] , (21) 

where S is the total spin of the skyrmion with respect to the uniform polarized background. In the case of the vortex 
S = S vort ex = l/2(2n + 1). 

4 Dual formulation of the model 

To get further insight into statistical interactions of vortices we will perform Hubbard-Stratanovich transformation on the 
model (^J). The phase gradients squared in the Lagrangian (^]) can be rewritten as 

4expz / Sx [-£± dkX Ad kX A] = I [Dl£]ezpi / fx - , (22) 



2m J J 2mpA 

where the auxiliary fields l£ have been introduced. In the next step the phases can be split into two parts xa = Xa^^a- X°' s 
are multivalued phases due to vortices/skyrmions while 77's are singlevalued components. At the present stage the Lagrangian 
is linearized in phase gradients. Functional integration over 77's leads to the conservation laws 

p A + O k [ I * +PA{a t + eAk) }=0. (23) 
m 



These conservation laws will be identically satisfied when we introduce a pair of dual gauge fields 

{ pa , + PA(a£ + eA k )] } = e^d u B? ^ \eT-H^ . (24) 



rn 



This definition together with the definition of the topological currents, 2ttK^ — £^ va d v d a XA-, leads, after some rearrangement 
and integration by parts, to the following dual Lagrangian 

L d=L b + ^e^a^a. - Yy^^ + e \) d » B * + ^K K a\ ■ (25) 

A 

Lb is the Lagrangian of the dual gauge field 

Let us concentrate now on the minimal coupling between the topological current and the dual gauge field. The topological 
current due to the antivortex/antiskyrmion moving along the trajectory is 

K° A = -S lA S 2 [x- £(*)] , K\ = -S lA i k (t)6 2 [x-at)} . (27) 



Only the topological current of the J, component is nonzero. It is like a current of a negatively charged point particle. What 
is the dual gauge potential which couples to the topological current? We know that the relation between the dual magnetic 
field and the original density is pa = —H^. Far from any topological defects the dual magnetic field is just H^i = — $lAPo 
and is directed down the z-axis. A single vortex/skyrmion moves in a uniform dual magnetic field. This is the origin of the 
Magnus force. Thus Magnus force is in fact a dual Lorentz force. Topological defects distort the uniform background. Their 
contribution to the dual gauge potential is 



where v runs over topological defects. With this form of the gauge potential we are able to work out the mutual statistical 
interaction 

-2n [ <Px K^Bjt = ^ Wd 2 . e M t5^{x-i w ) = -9k £ *Arg[& - Q (29) 

which is the same as Eq.(^ij). 
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Even if the inertial mass of quasiparticles is zero there still remains some possibility of motion provided that there are 
long range potential interactions between them. For nonzero Zeeman coupling the quasiparticles are exponentially localized 
field configurations. There is thus no long-range potential interaction through the matter fields. However there is still the 
long range Coulomb interaction (|^) which gives rise to the following effective Lagrangian for diluted quasiholes 

V v 1 v<w I St) Sic I v <w v 1 

The last term is a total time derivative and can be neglected in classical considerations. According to the classical part of 
this Lagrangian the particles move along the trajectories of constant total potential energy. For example a pair of quasiholes 
performs a strictly circular clockwise motion. 



5 Remarks 

We have derived statistical interaction of skyrmions in the framework of the time-dependent Ginzburg-Landau model. These 
interactions can be much stronger than interactions of fully polarized vortices. This conclusion could not be obtained in 
the effective Heisenberg model. This example shows that the Heisenberg model is not relevant not only to description of 
quasiparticles' energies but it also fails to give a correct answer to more basic topological questions. 

We did not make any quantitative predictions as to what is the skyrmion's spin S which determines the strenght of the 
statistical interaction. Such calculations are possible in the Ginzburg-Landau model but, there exist already microscopic 
Hartree-Fock techniques to establish the value of S, see . 

From our calculations it turns out that that the geometrical phase picked up by electronic wave- function during a clockwise 
exchange of two skyrmions is 

r -^ + 2(2^^^ + (^Tl)] , (3D 

where TVf is the number of electrons with reversed spins trapped inside the skyrmion core. This phase determines the quantum 
statistics of skyrmions. There are two types of anyons. For iVf even the phase is ^n+i) U P ^° an even niultipilicity of 7r. This 
is the same kind of statistics as for fully polarized Laughlin quasiparticles (or quasiholes). Condensation of such skyrmions 
leads to hierarchy states with odd denominator filling fractions. For skyrmions with iVj odd the exchange phase ( |3l| ) is 
^fn+f U P ^° unrelevant even multiplicities of it. This type of anyons differs from the Lauglin quasiparticles by an addition 
of one flux quantum. Their condensation gives rise to even-denominator filling fractions. To summarize the condesation of 
skyrmions of the primary state with v = 2rt 1 , 1 gives rise to the following filling fractions labelled by p 

p= 1,2,3.... , (32) 



(2n + 1) - 



where a = +1 or — 1 dependent on whether skyrmions or antiskyrmions condense. Note that in the above formula p is an 
arbitrary integer so that the filling fractions can have both even and odd denominators. 
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